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Abstract We present a calculation of the three-quark 
core contribution to nucleon and Z\-baryon masses and 
A electromagnetic form factors in a Poincare-covariant 
Faddeev approach. A consistent setup for the dressed- 
quark propagator, the quark-quark, quark-'diquark' and 
quark-photon interactions is employed, where all ingre- 
dients are solutions of their respective Dyson-Schwinger 
or Bethe-Salpeter equations in a rainbow-ladder trun- 
cation. The resulting A electromagnetic form factors 
concur with present experimental and lattice data. 

Keywords Nucleon • Delta • Form factors • Dyson- 
Schwinger equations • Bethe-Salpeter equation 

1 Introduction 

The exploration of the rich structure of the nucleon rep- 
resents one of the main tasks of contemporary particle 
physics. Present experimental facilities report accurate 
measurements of the nucleon's electromagnetic form 
factors. The lowest-lying excited state of the nucleon, 
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the Z\(1232) baryon, is produced at energies above the 
pion-production threshold and plays an important role 
in nuclear strong interactions. A comprehensive study 
of Z\-baryon properties in connection to those of the 
nucleon is expected to answer important issues of con- 
temporary research such as the chiral cloud content of 
baryons and its impact on baryon properties. In this 
view, a comparative analysis of the N A'y and the AA^ 
transitions will elucidate the nature of the Z\-baryon 
as a pure quark state, rather than a molecular state, 
and reveal connections between experimental observa- 
tions and the fundamental phenomena that govern the 
physics of hadronic constituents. 

A quark-core analysis of nucleon and A masses and 
electromagnetic form factors has recently been carried 
out in the Dyson-Schwinger approach [T]-[5]. Dynamical 
chiral symmetry breaking and confinement, two gen- 
uinely non-perturbative phenomena tightly connected 
with the formation of bound states, can be consistently 
addressed only within a non-perturbative approach to 
QCD. Such a framework is provided by the Dyson- 
Schwinger equations (DSEs) which are an infinite set 
of coupled integral equations for QCD's Green func- 
tions; see e.g. [BHS] for reviews. In this context, hadron 
properties are studied via covariant bound-state equa- 
tions [5HTT] : mesons {qq bound states) can be described 
by solutions of their Bethe-Salpeter equations (BSEs); 
baryons {qqq bound states) are studied by means of a 
covariant Faddeev equation. They are homogeneous in- 
tegral equations for a hadron's amplitude and depend 
on the dressed quark propagator as well as the quark- 
antiquark or three-quark kernel, respectively. 

The covariant Faddeev equation was recently solved 
for the nucleon mass by implementing a rainbow-ladder 
(RL) truncation, i.e. a dressed gluon-ladder exchange 
kernel between any two quarks [T^Hn| . While this puts 
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the analysis of baryon properties on the same footing as 
sophisticated meson studies, the numerical efforts are 
involved. To simplify the problem, a quark-'diquark' 
bound-state BSE to study baryon properties has of- 
ten been employed, see e.g. [T5]. It is based on the ob- 
servation that the attractive nature of quark-antiquark 
correlations in a color-singlet meson is also attractive 
for 3c quark-quark correlations within a color-singlet 
baryon ^^^7\. 

In connection with hadronic bound-state equations, 
the RL truncation of DSEs has been widely employed 
for studying hadron observables. It provides a reason- 
able description of pseudoscalar-meson, vector-meson 
and nucleon ground-state masses and electromagnetic 
properties, see e.g. [H [TSH21]. Other quantities, most 
notably the masses of axial-vector and pseudoscalar 
isosinglet mesons, are not reproduced so well. Efforts 
to go beyond RL have been made (see e.g. j2M27) ) but 
typically require a significant amplification of numeri- 
cal effort. In recent studies a consistent implementation 
of additional structures in the quark-gluon vertex and 
quark-antiquark kernel has proven capable to provide a 
better description of such observables as well pSH5D] . 

On the other hand, substantial attractive contribu- 
tions beyond RL come from a pseudoscalar meson cloud 
which augments the 'quark core' of dynamically gener- 
ated hadron observables in the chiral regime, whereas 
it vanishes with increasing current-quark mass. A view- 
point explored in Ref. [3T] was to identify RL with 
the quark core of chiral effective field theory which, 
among other corrections, must be subsequently dressed 
by pion-cloud effects. Such a quark core can be modeled 
by a current-quark-mass dependent input scale which 
is deliberately inflated close to the chiral limit. Re- 
sulting mass-dimensionful t:, p, N and A observables 
were shown to be consistently overestimated and mostly 
compatible with quark-core estimates from quark mod- 
els and chiral perturbation theory [21 HI [3T] , a pattern 
also present in a recent exploratory study of the QCD 
chiral transition temperature in this approach [32j . 

In the present work we follow this point of view 
and report on the latest calculations of the quark-core 
contributions to the A mass and electromagnetic form 
factors. 

2 Quark-diquark Faddeev-equation framework 

In a quark-'diquark' scenario, a color-singiet baryon 
emerges as a bound state of a color-triplet quark and 
color-antitriplet diquarks. Diquark correlations are im- 
plemented as a separable sum of pseudoparticle pole 
contributions in the quark-quark scattering matrix. This 
procedure simplifies the covariant Faddeev equation to 




Fig. 1 (Color online) The quark-diquark BSE. 



a quark-diquark BSE on the baryon's mass shell which 
is diagrammatically represented in Fig. [l] 

The baryon mass and amplitude are obtained as nu- 
merical solutions of the quark-diquark BSE once all its 
ingredients are specified: the dressed-quark propagator 
(single line), the diquark propagator (double line), and 
the diquark amplitude and its charge-conjugate P'^ 
which appear in the quark-diquark kernel. The binding 
mechanism in the baryon is realized via an iterated ex- 
change of roles between the single quark and any of the 
quarks contained in the diquark. 

The fundamental building block which connects the 
properties of baryons with the underlying structure of 
QCD is the renormalized dressed quark propagator. It 
involves the quark mass function MijP') which is non- 
perturbatively enhanced at small momenta and thereby 
indicates the dynamical generation of a large constituent- 
quark mass. This manifestation of dynamical chiral sym- 
metry breaking emerges in the solution of the quark 
DSE. The latter involves the gluon propagator and the 
quark-gluon vertex that both satisfy their own DSEs 
which depend on higher-order Green functions. In prac- 
tical calculations, the resulting infinite set of coupled 
DSEs is circumvented by employing truncations that 
preserve the underlying symmetries of QCD. 

In connection with meson properties, e.g. to estab- 
lish the pion as the Goldstone boson of spontaneous chi- 
ral symmetry breaking, it is imperative to preserve the 
axial-vector Ward-Takahashi identity. It connects the 
kernel of the quark DSE with that of the pseudoscalar 
meson BSE, ensures a massless pion in the chiral limit 
and leads to a generalized Gell-Mann-Oakes-Renner 
relation [331 IM]- The simplest truncation that satis- 
fies this constraint is the rainbow-ladder (RL) trunca- 
tion which amounts to an iterated dressed-gluon ex- 
change between quark and antiquark and has been ex- 
tensively used in Dyson-Schwinger studies of hadrons, 
see e.g. [22l [24] and references therein. 

The RL truncation retains only the vector part of 
the dressed quark-gluon vertex. Its non-perturbative 
dressing, together with that of the gluon propagator, 
is absorbed into an effective coupling a{k'^) which rep- 
resents the only unknown function of the model. Herein 
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we employ for a(fc^) the frequently used ansatz [35] 

a{k') = n^' (^^^e-^l^) + aMk') , (1) 

where k is the gluon momentum. At large gluon mo- 
menta, the second term in the effective coupling a(fc^) 
decreases logarithmically and reproduces QCD's per- 
turbative running coupling. At small and intermedi- 
ate gluon momenta, the first term must exhibit suffi- 
cient strength to allow for dynamical chiral symmetry 
breaking and the dynamical generation of a constituent- 
quark mass scale. The infrared behavior of the effective 
coupling is controlled by two parameters: an infrared 
scale A and a dimensionless width parameter rj. 

The features of the effective coupling directly trans- 
late from the case of mesons to that of diquarks which 
enter the quark-diquark bound-state equation in Fig.[l] 
By virtue of the RL truncation, both mesons and di- 
quarks are bound by the same gluon-exchange mecha- 
nism. The corresponding diquark BSEs are obtained by 
assuming timelike diquark poles at certain values of the 
total diquark momentum in the quark-quark scat- 
tering matrix, i.e. = —m^^, P^ = —m'^^, which char- 
acterize the lightest diquarks, namely the scalar and 
axial-vector ones. Diquarks carry color and are hence 
not observable; yet such a pole structure does not con- 
tradict diquark confinement, see e.g. [7]. In the present 
context, timelike diquark poles emerge as an artifact 
of the RL truncation which does not persist beyond 
RL [53]; nevertheless they indicate the presence of di- 
quark mass scales within a baryon. Studies in a similar 
setup provide further support for diquark correlations 
as a reasonable concept for the description of baryons 
[37] • While both scalar and axial- vector diquark corre- 
lations are important for the description of the nucleon, 
the spin— 3/2 and isospin— 3/2 flavor symmetric A ne- 
cessitates only axial- vector diquark correlations. 

The scalar and axial-vector diquark BSEs only spec- 
ify the on-shell diquark amplitudes whereas diquarks 
in a baryon are off-shell. Within the separable ansatz 
for the scattering matrix, information on its off-shell 
behavior is encoded in the scalar and axial-vector di- 
quark propagators. By reinserting the separable pole 
ansatz into the Dyson series for the scattering matrix, 
the resulting diquark propagators are completely spec- 
ified from their substructure; see [2j [11] for details. 

3 Nucleon and A masses 

All the ingredients of the quark-diquark BSE are now 
determined: the quark propagator is obtained as a so- 
lution of the quark DSE, the scalar and axial-vector 



diquark amplitudes as solutions of the diquark BSEs, 
and the respective diquark propagators follow from the 
separable diquark-pole ansatz. All these elements enter 
the quark-diquark kernel in Fig. [T] and are numerically 
calculated within RL truncation which involves only 
one parametrization as its input, namely the effective 
coupling a(/c^) of Eq. ([T]). 

Upon a decomposition of the quark-diquark ampli- 
tudes of nucleon and A into orthogonal sets of Dirac 
covariants, their amplitudes and masses emerge as nu- 
merical solutions of the respective quark-diquark BSEs 
(details on the calculation were reported in Ref. [2]). 
We depict the results for nucleon and A masses, to- 
gether with that of the p-meson, in Fig. [2] and com- 
pare their evolution with to lattice calculations. 
The masses are calculated within two versions of the 
model which are characterized by the infrared scale A 
in the effective coupling of Eq. ([T]). The dashed lines 
in Fig. [2] are obtained by using a fixed scale A = 0.72 
GeV which is adjusted to reproduce the experimental 
pion decay constant and kept fixed for all values of the 
quark mass. It yields the result = 0.94 GeV and 
Af/i = 1.28 GeV at the physical u/d-quark mass corre- 
sponding to a pion mass = 140 MeV. These results 
are reasonably close to the experimental values and 
consistent with pseudoscalar-meson and vector-meson 
ground-state properties which are satisfactorily repro- 
duced in this setup, e.g. [2l|20j|24], and moreover insen- 
sitive to the shape of the coupling in the infrared, i.e. 
to a variation of the width parameter rj |35) . 

In the second version of the model, the hadronic 
quark-core properties are implemented through a current- 
quark mass dependent scale A in Eq. ([T]) . It is deliber- 
ately inflated close to the chiral limit and flxed to repro- 
duce the core properties of the p meson 31j . The corre- 
sponding results in Fig. [2] are depicted by bands which 
indicate the variation with the width parameter rj. For 
a value of yl = 0.98 GeV at the u/d mass, the resulting 
values are M^°'''= = 1.26(2) GeV and M^°'''= = 1.73(5) 
GeV. In the chiral region, the core version of the model 
uniformly overestimates the experimental and lattice 
data as well as the results obtained using the fixed- 
scale model, whereas this deviation decreases with in- 
creasing current-quark mass. While pionic corrections 
to hadronic observables vanish at large quark masses, 
they should reduce the core masses of nucleon and A by 
~ 300 MeV in the chiral region. In this respect, the core 
value for is roughly consistent with a pseudoscalar- 
meson dressing providing the dominant correction to 
the quark-diquark core, whereas the somewhat large 
deviation between the experimental and 'core' mass of 
the A may indicate the relevance of further diquark 
channels in describing A properties. 
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Fig. 2 Evolution of p-meson, nucleon and A masses with 
as obtained from their qq and quark-diquark BSEs. The bands 
correspond to the core model and represent the sensitivity of the 
masses to the width parameter r] = 1.8 it 0.2. The dash-dotted 
lines depict the results in the fixed-scale model for the central 
value r] = 1.8. Stars denote the experimental values |38| of p, N 
and A at the u/d— quark mass and 0, i7 at the strange-quark mass 
whose positions are indicated by vertical lines. Note that there is 
no ss pseudoscalar meson in nature; the value m^s = 0.69 GeV 
corresponds to a meson-BSE solution at a strange-quark mass 
tUb = 150 MeV |34j . We compare with a selection of lattice data 
for mp [39] Ho], Mjv ]4Tti44] and [iilliS] . 



We note that simple relations between the two se- 
tups hold in the chiral limit [TT] where the infrared pa- 
rameter A in Eq. ([T]) represents the only relevant scale 
in the system, i.e. the scale of dynamical symmetry 
breaking. This implies that mass-dimensionful quanti- 
ties which are sensitive to the infrared properties scale 
with A. For a given set of dimcnsionful observables, 
the fixed-scale model produces results that are overes- 
timated by the same percentage upon entering its core 
version. On the other hand, the distinction between di- 
mensionless quantities calculated within either versions 
of the model becomes irrelevant in the chiral region. As 
we will point out below, this is the case for electromag- 
netic form factors. 



4 Delta electromagnetic form factors 

To compute the electromagnetic form factors of the A- 
baryon, one must specify how the photon couples to its 
constituents. In the quark-diquark context this amounts 
to resolving the coupling of the photon to the dressed 
quark and the axial-vector diquark (impulse approxi- 
mation), to the axial- vector diquark amplitude (seag- 
ulls) and to the exchanged quark in the quark-diquark 
kernel. With this decomposition the A electromagnetic 
current is automatically conserved [46] . 

At the level of the constituents, current conservation 
translates to electromagnetic Ward-Takahashi identi- 
ties (WTIs) which constrain the longitudinal parts of 
the above vertices and unambiguously relate them to 
the previously determined quark and diquark propaga- 
tors and diquark amplitudes. On the other hand, cur- 
rent conservation only partly constrains pieces which 
are transverse to the photon momentum. Their impor- 
tant role in physical observables is accounted for by 
augmenting the vertices as determined from their WTIs 
by appropriate transverse /o-meson pole contributions. 
The details of the construction of the electromagnetic 
current are presented in Appendix C of [1]. 

Following these prescriptions for the electromagnetic 
current, the A electromagnetic form factors — the elec- 
tric monopole Geo{Q^), electric quadrupole G_e2(Q^), 
magnetic dipole G'mi(Q^) and magnetic octupole form 
factor Ga[3{Q^) — are directly related to the effective 
quark-gluon coupling in Eq. Q. In Fig. |3] we depict 
the 'core' contributions to the A'^ electromagnetic form 
factors and compare to recent lattice data [H]. Due to 
isospin symmetry the A'^'^ , A'^ and A~ form factors are 
simply obtained by multiplying those of the Z\+ with 
the appropriate charges. The electromagnetic form fac- 
tors (as dimensionless quantities) are plotted as a func- 
tion of the dimensionless variable Q^/M^ for our data, 
and / (Af^*)^ for the lattice data. This enables an un- 
ambiguous comparison between our form factor results, 
calculated with an implicit 'core' mass Ma > M'^^ and 
the lattice results. From another perspective, the cal- 
culated Ma sets the scale of dynamical chiral symme- 
try breaking in either version of the model, and hence 
the dimensionless form factors in Fig. [3] would approxi- 
mately match the corresponding dimensionless quanti- 
ties calculated within the fixed-scale model. 

Experimentally only the magnetic moments of the 
A^ and Z\^~'' are known, albeit with large errors. For 
instance, for the A'^ the Particle Data Group quotes 
the value S.S^l^^ g [3H]. Our result for the magnetic mo- 
ment, GAfi(O) = 3.64(16), compares well with quark- 
model predictions and chirally extrapolated lattice re- 
sults; see [37] and references therein. The deformation 



5 




0.5 1.0 



1 .5 0.0 




0.5 1.0 





\MeV\ 




691 




509 


\ 


384 




1.5 0.0 




Fig. 3 (Color online) Electromagnetic form factors of the A calculated at the physical point T7i,r = 140 MeV, and compared to 
unquenched lattice data of Ref. |44| at three different pion masses. The bands represent the sensitivity to a variation of the infrared 
width parameter = 1.8 it 0.2. Adapted from Ref. 1 



of the Z\-baryon is encoded in its electric quadruple 
and magnetic octupole moments. Currently there are 
no experimental measurements for these observables. 
From the measurement of the N^A transition, one can 
extract the value Ge2{0) = —1.87(8) in the large- A^c 
limit [milH]; similar values are predicted by constituent- 
quark models [49] . Lattice calculations indicate a nega- 
tive value for Ge2{0) as well j44] but are limited by large 
statistical errors. Our result for the electric quadrupole 
moment, Ge2{0) = —1-32(16), is negative and com- 
patible with these observations. We note that Ge2{Q^) 
develops a zero-crossing at Q'^/M'^ ~ 0.6, a feature 
which is unexpected but not clearly excluded from the 
available lattice results. Our calculation for the mag- 
netic octupole moment yields a small and negative value 
Gm3{0) = —0.26(4). We note that the electric quadrupole 
and magnetic octupole form factors are negative through- 
out the current-quark mass range which indicates an 
oblate deformation of the Z\'s charge and magnetiza- 
tion distributions. 

We conclude that the rainbow-ladder truncated Poin- 
care-covariant Dyson-Schwinger/Bethe-Salpeter setup. 



upon implementing an appropriate input scale, pro- 
duces consistently overestimated core contributions to 
nucleon and A masses. The overall attractive effect of 
chiral corrections is expected to shift the core masses to 
the experimental values. A 20 — 30% reduction for dy- 
namically generated hadron masses in the chiral region 
is anticipated, while pionic effects decrease at larger 
quark masses. 

The impact of chiral corrections upon the low-Q^ 
behavior of the A electromagnetic multipole form fac- 
tors remains at present unclear. While near-future mea- 
surements at MAMI and JLab facilities remain to vali- 
date our predictions for the Z\'s electromagnetic prop- 
erties, the results collected herein show good agreement 
with lattice results and quark model analyses. Towards 
a complete understanding of the structure of baryons, 
our approach can be improved by implementing chiral 
corrections to the rainbow-ladder truncation, and aug- 
mented by studies which eliminate the diquark ansatz 
in support of a fully Poincare-covariant solution of the 
three-quark Faddeev equation. Fruitful insight can be 
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gained by a forthcoming extension of our approach to 
the investigation of the NAj transition. 
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